Density-dependent deformed relativistic Hartree-Bogoliubov theory in continuum 
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The deformed relativistic Hartree-Bogoliubov theory in continuum with the density-dependent 
meson-nucleon couplings is developed. The formulism is briefly presented with the emphasis on 
handling the density-dependent couplings, meson fields, and potentials in axially deformed system 
with partial wave method. Taking the neutron-rich nucleus 38 Mg as an example, the newly developed 
code is verified by the spherical relativistic continuum Hartree-Bogoliubov calculations, where only 
the spherical components of the densities are considered. When the deformation is included self- 
consistently, it is shown that the spherical components of density-dependent coupling strengths are 
dominant, while the contributions from low-order deformed components are not negligible. 
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Since the experimental discovery of a large neutron 
radius in 11 Li the exotic nuclear halo phenomenon 
becomes one of the most interesting topics close to the 
nucleon drip lines. In order to describe the halo phe- 
nomenon, the asymptotic behavior of nuclear densities 
at large distance from the center must be treated prop- 
erly, and the discrete bound states, the continuum, and 
the coupling between them need to be dealt with simul- 
taneously in a self-consistent way. 

During the past decades, the covariant density func- 
tional theory (CDFT) has achieved great success in de- 
scribing lots of nuclear phenomena in both stable and ex- 
otic nuclei , including the recent achievements in nu- 
clear magnetic moments 1(1 , pseudospin symmetry @, [|| , 
low- lying excitations [9l4llj. magnetic and antimagnetic 
rotation fl2T4l4l ]. collective vibration [l5l - [l9j . and so on. 

In particular, great efforts have been dedicated 
to dev eloping the relativistic Hartree-Bogoliubov 
(RHB) [20-24J and relativistic Hartree-Fock-Bogoliubov 
(RHFB) [25| theories in continuum for a self-consistent 
description of spherical halo nuclei. In order to describe 
the halo phenomena in deformed nuclei, a deformed RHB 
theory in continuum has been developed recently [26l — 
[30l | . An interesting shape decoupling between the core 
and halo in 42 ' 44 Mg has been found [H, [29|. In these 
applications, the deformed RHB equations are solved in 
a Woods-Saxon basis 
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with the partial wave method, 
and only nonlinear meson self-coupling interactions are 
used so far. 

In recent years, the RHB models with the density- 
dependent meson-nucleon couplings have attracted more 
and more attention owing to improved descriptions of 
the equation of state at high density, the asymmetric nu- 
clear matter, and the isovector properties of nuclei far 



from stability [32M37l |. Therefore, it's necessary to de- 
velop a density-dependent deformed relativistic Hartree- 
Bogoliubov (DDDRHB) theory in continuum. 

In this brief report, the deformed RHB theory in con- 
tinuum with density-dependent meson-nucleon couplings 
is developed for a wider compatibility of modern func- 
tionals. The key technique here is to handle the density- 
dependent couplings and the meson fields as well as the 
potentials in a deformed system with the partial wave 
method. 

The starting point of the CDFT is a Lagrangian den- 
sity where nucleons as Dirac spinors interact with each 
other by exchanging effective mesons (a, uj, and p) and 
photons 0-1U, 
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where M is the nucleon mass, and m a ,g a , m u , g u , m p , g p 
are the meson masses and density-dependent coupling 
strengths of the respective mesons. The field tensors for 
the vector mesons and photons are 



R» v 
F^ 



dPA" -FA*. 
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Following the formalism in Ref. [38(, one can derive the 
RHB equation, 
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with the quasiparticle energy E/-, Fermi surface A, and 
Dirac Hamiltonian 

h D = a • p + V(t) + 0[M + S(t)] + S rca r(r), (4) 

where the Dirac matrices a = 7°7 and /3 = 7 . The 
scalar and vector potentials are respectively, 



with 



S(r) = .g ff cr(r), 



(5a) 



V(r) = g u 0J (r) + g p T 3 po(v) + e (1 ^ A (v), (5b) 
and the rearrangement term is 

Srear(r) = ^ (r)<r(r) + ^ Pv (v)m (r)+^p 3 (r)po « 
op v op v dp v 

(6) 

which comes from the densit y-dep endent behaviors of the 
meson- nucleon couplings [35l . |36[ , 

[ 9V>(p sa t) exp[-a0(a; - 1)] for = p, 

where p sat denotes the baryonic saturation density of nu- 
clear matter and x = p v j Psat • 

The equations of motion for mesons and photons read 



(-A + m )</>(r) = s (r), 



(8) 



where to^ are the meson masses for <p = cr, w, p and zero 
for the photons. The corresponding source terms are 



s 4 r ) = 



-9a(Pv)Ps{r), for the a field, 

9u(Pv)Pv(t), for the u field, 

9 P {Pv)p3(r), for the p field, 

ep c (r), for the Coulomb field 
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with various densities 



Ps(r) 


k>Q 


(10a) 


Pv(r) 


= ]T^(r)^(r), 

fc>0 


(10b) 


P3(r) 


= ^^(r)r 3 ^(r), 

fc>0 


(10c) 


Pc(r) 


= E^w^^w- 


(lOd) 



fe>0 

calculated within the no-sea approximation. The sum 
over k > runs over the quasiparticle states correspond- 
ing to single-particle energies in and above the Fermi sea. 

For axially deformed nuclei with spatial reflection 
symmetry, the potentials in Eqs. ([5]) and (j6)), coupling 
strengths in Eq. ([7}, meson fields in Eq. ©, and densi- 
ties in Eqs. (1101) are expanded in terms of the Legendre 
polynomials |39jj, 



f(r) = ^2f x (r)P x (cos6), A -0,2,4, 



(11) 
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d(cos6»)/(r)P A (cos6>). (12) 



According to the partial wave method, first, the den- 
sities in Eqs. (fTO)) are represented as 



(13) 



(14) 



p(r) = Px(r)P\ (cos 6). 

X 

Second, the coupling strengths in Eq. ([7]) read 

94>{Pv) = J29<l>A r ) p x(cos9). 
x 

By taking the a meson as an example, 
2A + 1 f 1 

9<j,\{r)= — - — J d(cos9)P\(cos9)g r7 (p sat ) 

1 + MKai Pv,\i{r)P\i (cos 9)}/p S3 ,t + dy} 2 

CT 1 + c <t{Ea2 /°f 1 A2(r)PA2(cOS0)]//O 8a t + C^} 2 ' 

(15) 

as illustrated for the neutron-rich nucleus 38 Mg in Fig. [3] 
below. Third, the meson and Coulomb fields are, re- 
spectively, solved by the Klein- Gordon and Poisson equa- 
tions 15||, i.e., 



(r) = J dv'D{r,e,r\e'-m^)s^{r\e'), 

= Y,Mr)Px(cos0), (16) 



with 

Mr) 



hxiim^r) I dr' V' 2 j\{im^r')s<i,^(r') 
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+ jxiim^r) I dr'r' hxiim^r^s^^r') 



for the meson fields, where j\ and h\ are the spherical 
Bessel and Hankel functions in static Green functions 
D(r, 8, r', 9'; m^), and 



1 r r°° l 



(18) 



for the Coulomb field. Finally, the scalar potential S(r) 
in Eqs. ([5]) is written as 



S(r)=^S A (r)P A (cos< 



(19) 



with 
Sx(r) 



2A+ 1 



E 9a,\i{r)crx2(r) 



X1X2 



d(cos (9)P A i (cos 6)Px2 (cos 0)P a (cos 0), 



(20) 
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where the partial waves Ax and A2 arc coupled to A. 
The vector potential V(r) is obtained in the similar way. 
It should also be noted that the rearrangement term in 
Eq. ©, 



(r) = Z Icai ,\(r)P\(cos0), 



(21) 



is a result of the couplings among three partial waves. 
The total energy of a nucleus is composed of 



E u 



E 



part 



E a + E u + E p + E IC 



E CO ul + E c 



(22) 



where the center-of-mass correction E cm can be calcu- 
lated in an empirical or microscopic way [29l l36l |40| - |42| . 

In present Brief Report, we use the functional PKDD 
developed in Ref. [36| for the particle-hole channel and a 
zero range density-dependent pairing force, 



V pp {r, r') 



^ (1 _r)(l-%(r-r'), (23) 

<S Psat 



for the particle-particle channel. The pairing strength 
V = 348.40 McV fm 3 and energy cutoff E^ u p t = 60 
MeV in the quasiparticle space are fitted for reproducing 
the proton pairing energy in the spherical nucleus 20 Mg 
given by the RHB calculation with Gogny pairing force 

dis m. 

Following the corresponding convergence check in 
Ref. [1^], the DDDRHB equations are solved in a spher- 
ical Dirac Woods-Saxon basis [3l[ determined by a box 
size i? max = 20 fm and a mesh size A?' = 0.1 fm with 
the energy cutoff E^ ut = 150 MeV and quantum number 
cutoff for angular moment j m &x = 19/2. For each 
block in Woods-Saxon basis, the number of negative en- 
ergy states in the Dirac sea is the same as that of positive 
energy states. 

Furthermore, the cutoffs are also necessary for the Leg- 
endre expansions of the densities, coupling strengths, 
meson fields, and potentials in the calculations. They 
are denoted with A P!inax , A s , max , A^, max , and Xv,ma^, re- 
specti vely, w here A p , max = Xv,max = 4 arc used as in 
Refs. [28|, [29j. in the following, we investigate the depen- 
dence of the DDDRHB results on A g>max and A^ max due 
to the density-dependent meson-nucleon couplings. 

In Fig. [TJ the total energy E to t, matter root mean 
square (rms) radius R m , and quadrupole deformation (3 
are plotted as functions of the cutoffs in Lcgendre expan- 
sions for the coupling strengths Ag imax and meson fields 
A^max, respectively. The convergence of these quanti- 
ties with Ag jmax and A0 imax is seen. The relative differ- 
ences of -Etot, Rm, and {3 between the calculations with 



g.max 



10 (A , r 



6) and A 9j 
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A 

are less than 0.001%. Therefore, the cutoffs A g>max = 10 
and A0 imax = 6 are adopted in the following calculations. 

In order to verify the accuracy of the present DDDRHB 
code, we first calculate the neutron-rich nucleus 38 Mg 
with DDDRHB but constrained to the spherical case by 
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FIG. 1: (Color online) Total energy Etot, matter rms radius 
R m , and quadrupole deformation /3 of 3S Mg as functions of the 
cutoffs in Legendre expansions for (a) the coupling strengths 
A 9ima x and (b) the meson fields A^.max, respectively. 



TABLE I: Ground state properties of Mg calculated by 
RCHB [H| as well as DDDRHB in the spherical case 
(A Pl max = 0) and deformed case (A p , ma x = 4) with the func- 
tional PKDD [3(|. The rms radii for neutron R n , proton R p , 
and matter Rm are in units of fm. The total energy E to t 
and the corresponding contributions from the particles E pal t, 
meson fields E a , E u , and E p , Coulomb field -B cou i, rearrange- 
ment term -Erear, pairing -B P ai r , and center-of-mass correction 
E cm are in units of MeV. 



RCHB 



DDDRHB 







Agjmax — 


hp, max — 4 


p 






0.46 


Rm 


3.73 


3.73 


3.80 


Rp 


3.04 


3.04 


3.12 


Rm 


3.53 


3.53 


3.60 


-Epart 


-856.36 


-856.32 


-877.05 


Eu 


4886.42 


4887.47 


4919.31 


E u 


-4070.70 


-4071.61 


-4103.55 


E P 


-38.70 


-38.73 


-38.90 


Ecoul 


-32.20 


-32.20 


-32.02 


Ei-car 


-121.29 


-121.30 


-119.45 


Epair 


-14.53 


-14.57 


-1.67 


Ecm 


-7.25 


-7.25 


-7.80 


Etot 


-254.61 


-254.50 


-261.14 



taking A Pimax = 0. The obtained bulk properties includ- 
ing the rms radii as well as the total energy and con- 
tributions from each component are listed in the middle 
column of Table HI For comparison, the corresponding 
results calculated by the spherical relativistic continuum 
Hartree-Bogoliubov (RCHB) [23| theory are shown in the 
left column. It is found that both calculations agree with 
each other quite well. The rms radii are the same up 
to 0.01 fm, and the total energies are the same up to 
0.11 MeV which corresponds to an accuracy of ~ 0.04%. 
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FIG. 2: (Color online) Single-particle energies E ca , n of neu- 
trons in the canonical basis above —20 MeV and their occu- 
pation probabilities v 2 of 38 Mg calculated by (a) RCHB and 
by DDDRHB in the (b) spherical case (A p , ma x = 0) and (c) 
deformed case (A p , ma x = 4). The Fermi surfaces are shown as 
the dotted lines. 



In addition, the neutron single-particle energies in the 
canonical basis of 38 Mg calculated by DDDRHB with 



A^max = are compared to those by RCHB in Fig. [5] 
The length of each level is proportional to the occupation 
probability v 2 . It can be seen that an excellent agreement 
is achieved. 

In order to investigate the deformation effect, the 
neutron-rich nucleus 38 Mg is calculated with DDDRHB 
by taking A Pimax = 4 following the convergence study in 
Ref. [26|. The corresponding bulk and neutron single- 
particle properties are also shown in Table U and Fig. [5J 
respectively. In this case, the minimum in the poten- 
tial energy surface locates at a prolate deformation with 
j3 = 0.46, and the rms radii increase by ~ 2%. Due to 
the spherical symmetry breaking, the degenerate single- 
particle levels are split and labeled by the third compo- 
nent of the total angular moment and parity f2 x . 

In order to illustrate the partial waves of density- 
dependent coupling strengths, we show in Fig. [3] the g a ,\ 
in Eq. (fT4")) for 38 Mg. It can be seen that the major com- 
ponent is that with A = 0, which is more than one order of 
magnitude larger than the others. The amplitudes of g a ^\ 
decrease quickly with increasing A and become negligible 
when A ^ 8. This confirms that the cutoff A g . max = 10 is 
reliable. 

In summary, the density-dependent deformed relativis- 
tic Hartrcc-Bogoliubov theory in continuum is developed. 
The key formalism on handling the densities, coupling 
strengths, meson fields, and potentials in a deformed 
system with the partial wave method is presented. The 



newly developed DDDRHB code is verified by compar- 
ing the bulk and single-particle properties of the neutron- 
rich nucleus 38 Mg calculated in the spherical case with 



A 



p,max 



to those obtained by the spherical RCHB. 



38 Mg is also studied by 
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corresponding minimum in the potential energy surface 
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As an illustration, the nucleus 
DDDRHB in the deformed case with A P;max = 4. 
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FIG. 3: (Color online) Decompositions of coupling strength 
g a .x with (a) A = and (b) A = 2, 4, 6, 8, 10 for 38 Mg. 



locates at a prolate deformation with j3 = 0.46. In this 
case, the major components of density-dependent cou- 
pling strengths are those with A = 0, while the contribu- 
tions from A = 2, 4, 6 components arc not negligible. 
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